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The concept of swim pressure quantifies the average force exerted by microswimmers on confining 
walls in non-equilibrium. Here we explore how the swim pressure depends on the wall curvature 
and on the presence of sharp corners in the wall. For active Brownian particles at high dilution, we 
present a coherent framework which describes the force and torque on passive particles of arbitrary 
shape, in the limit of large particles compared to the persistence length of the swimmer trajectories. 
The resulting forces can be used to derive, for example, the activity-induced depletion interaction 
between two disks, as well as to optimize the shape of a tracer particle for high swimming velocity. 
Our predictions are verifiable in experiments on passive obstacles exposed to a bath of bacteria or 
artificial microswimmers. 


There is a rapidly growing interest in understanding 
the physical principles of microswimmers both for liv¬ 
ing systems such as bacteria lljja and for artificial self- 
propelled colloidal particles [6l-l8|. Swimming of these 
active particles occurs at low Reynolds numbers and is a 
nonequilibrium phenomenon. Intriguingly, the concepts 
of equilibrium statistical physics break down for this new 
type of ” active matter”. 

In particular, when passive obstacles are exposed to 
an active suspension, such as a bacterial bath, they ex¬ 
perience non-thermal fluctuations from the self-propelled 
particles which result in a fascinating wealth of new 
nonequilibrium phenomena ranging from self-starting ro¬ 
tors and cogwheels |9|;ll| and the spontaneous motion 
of microwedges [II, [l3| to the rectification and sorting 
of bacterial motion through asymmetric barriers 14l4lfil| 
and ratchets Q. This opens the possibility to power mi¬ 
croengines by an active bath Elis and to steer micro¬ 
robots by activity [l^. Moreover, more recently, the 
swimmer-induced interaction between two parallel walls 
H-ili and two spheres [ 2 ^ has been considered and 
simulated. This is in a certain analogy to what is known 
as depletion and fluctuation-induced [ 2 ^ forces in 
equilibrium (or active i) suspensions. 

There is a need to understand the principles of these 
novel non-equilibrium effects within conceptually simple 
models which generalize the concept of equilibrium ther¬ 
modynamics. This is most directly done in the context 
of active Brownian motion describing self-propelled col¬ 
loidal spheres. One of the basic quantities characterizing 
the interaction between a swimmer bath and an obstacle 
is the recently introduced swim pressure [ 2 ^, [2^. The 
swim pressure is defined as the average force per area 
(or per length in two spatial dimensions) exerted by a 
microswimmer bath on a planar wall, such that an in¬ 
tensive mechanical quantity can be defined in nonequi¬ 
librium suspensions 28|, [29| . In the limit of planar walls 
and high dilution of spherical swimm ers, the swim pres¬ 
sure has been calculated analytically 26j, l27| . Addition¬ 


swimmers, as well as the behavior of swimmers in strong 


confinement 32 


ally, there has been recent interest in the forces exerted 
on curved walls by both spherical and elongated [3l| 


In this article, we consider the swim pressure on curved 
walls, including walls with sharp corners. In particular, 
we study different wall shapes exposed to dilute suspen¬ 
sions of torque-free, active Brownian particles. We ob¬ 
tain a coherent and unified framework to describe both 
the motion of obstacles of arbitrary hard shape in a bac¬ 
terial bath and the swimmer-induced force between two 
spheres in close proximity. Our superposition-like the¬ 
ory applies to any general shape with length scales much 
larger than the persistence length of the active particle 
trajectories. Note that the opposite limit, where swim¬ 
mers are trapped in confinements much smaller than their 
persistence length, was extensively studied by Fily et al 
[321 [33| . The theory presented here is semi-empirical, 
but captures the essential physics similar to the Asakura- 
Oosawa model for depletion forces between spheres in a 
thermal solvent (23 |. 

In two spatial dimensions, we confirm that the curva¬ 
ture correction to the swim pressure scales as 1/R for 
small curvatures and we propose a fit for this correc¬ 
tion for arbitrary curvature. This is in some formal anal¬ 
ogy to equilibrium thermodynamic for surfaces where the 
curvature corrections to the interfacial free energy are of 
prime interest and still debated di-lsii. The corners are 
treated as a singularity in the shape and studied sepa¬ 
rately. The total force is obtained from a superposition 
principle. As applications, we calculate the force acting 
on some hard particles of V-like and U-like shape, and 
propose a simple formula for the depletion force between 
two disks which generalizes the seminal Asakura-Oosawa 
model to an active suspension. Our predictions are ver¬ 
ifiable in experiments on passive obstacles of different 
shape exposed to a bath of bacteria or of artificial mi¬ 
croswimmers. 

In our model, we consider a two-dimensional system 
consisting of ideal (i.e. non-interacting) self-propelling 
particles, in contact with one or more flat or curved walls. 
Apart from the self-propulsion force /o, the particles ex¬ 
perience a drag force with friction coefficient C,, and ro- 
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tational diffusion of their axis of self-propulsion with a 
diffusion constant Dr- The equations of motions for the 
position X and orientation </> are therefore: 


X = iextK + fo/Cn = iext/C + Won 

(j) = y/2DrR{t), 


( 1 ) 


where /ext is any external force the particles experi¬ 
ence from e.g. the confining walls, vq = /o/C is the 
self-propulsion velocity, and i?(t) represents a delta- 
correlated stochastic noise term with zero mean and unit 
standard deviation. Note that the particles do not un¬ 
dergo translational (Brownian) diffusion. 

The only inherent length scale in this model is the per¬ 
sistence length of the particle trajectories Ip = uq /Dr', the 
typical distance a particle travels before it rotates sig¬ 
nificantly away from its original orientation. Additional 
length scales can be imposed by the shape of the walls in 
contact with the active bath. 

In several recent studies 2^ 13, 371, it was shown that 
for an ideal system of swimmers, the swim pressure on a 
flat wall in two dimensions is given by: 


Po = pC,vl/2Dr, 


( 2 ) 


where p is the bulk density of swimmers. This swim 
pressure has been shown to behave similar to a thermo¬ 
dynamic state function only if the swimmers experience 
no external torques and the swim speed is uniform and 
isotropic throughout the system [2^ . 

In order to investigate the effects of curvature on this 
swim pressure, we calculate the density profile of torque- 
free swimmers at and near walls that are either curved or 
contain a corner. To this end, we find the steady-state so¬ 
lution to the Smoluchowski equation for the orientation- 
dependent density distribution p(x, </>, t) in the system: 


dp{x,(j),t) 

dt 


-v(())) ■ Vp(x, 


w ■ 


Here, v(/)) = (uq cos </>, vq sin /p) is the propulsion velocity 
of a particle with orientation 4>. We find the steady- 
state solution of this equation numerically on a mesh, for 
various choices of the hard-wall geometry. In all cases, 
we use at least 500 spatial cells, and 64 discrete orien¬ 
tations for the particles. The cells are adapted to the 
expected density profile, with smaller cells near the hard 
boundaries and sharp corners, where rapid changes in lo¬ 
cal density occur. Note that at the walls, the local density 
technically diverges, as the swimmers have a finite prob¬ 
ability of being located at the wall. To circumvent this 
issue, the cells directly adjacent to the hard boundaries 
are chosen particularly small (thickness on the order of 
0.001/p). Particles in these cells are considered to be in 
contact with the wall, and therefore contribute to the 
pressure experienced by the boundary. We have verified 
that varying the number and sizes of cells does not signif¬ 
icantly affect our results. Note that to avoid non-physical 



FIG. 1. Pressure P on walls as a function of the wall radius 
of curvature R. The red points were obtained by numerically 
solving the Smoluchowski equation, the blue points from sim¬ 
ulations, and the black line is the fit from Eq. Note that 
different choices of Ip collapse onto the same curve. The insets 
show a typical particle trajectory for each geometry. 


flow of particles at the wall, the cell boundaries of these 
edge cells must be chosen normal to the wall. Further 
details are described in the Appendix. 

In addition to numerically solving the Smoluchowski 
equation, we also perform Brownian dynamics simula¬ 
tions of active solutions in contact with hard walls. In 
these simulations, the point-like active particles do not 
interact, and move according to the overdamped Brow¬ 
nian equations of motion (Eq. (HJ). The interaction be¬ 
tween the particles and the wall is a purely hard repul¬ 
sion. 

In order to investigate the force exerted on a curved 
wall by an active solvent, we study both active solutions 
in spherical confinement, and active solutions in contact 
with the outside of a circular object, for a range of values 
of the radius of curvature R of the hard boundary. As ex¬ 
pected, in the limit of large R^ Ip, the average pressure 
on the walls reduces to the swim pressure Pq in Eq. ©■ 
However, when the radius of curvature is finite, particles 
in contact with the inside of the boundary tend to spend 
a longer time in contact with the wall than those trapped 
at the outside. This results in a correction to the swim 
pressure which increases in magnitude as the radius of 
curvature decreases. This observation is consistent with 
recent work by Mallory et al., where the authors stud¬ 
ied the net motion of curved objects placed in an active 
bath [ 3 ^. In agreement with their results, we find that 
the first-order correction to the active pressure is propor¬ 
tional to 1/i?, with the same coefficient for positive and 
negative radius of curvature. In Fig. [TJ we show the force 
on both the inside (i? < 0) and outside [R > 0) of the 
wall, as obtained from both simulations and the numer¬ 
ical steady-state solutions of Eq. ®. We find excellent 
agreement between both methods. The resulting curve 
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is well fitted by the empirical expression; 


PjR) 

Po 


ae + (1 


X log 2 

log(l +exp[;p/i?])’ 


( 4 ) 


where a ~ 0.0893 is a single fitting parameter. Note 
that this expression reduces to 1 in the limit of flat walls 
(i? —oo), as expected. We emphasize here that the 
bulk swim pressure Pq is defined via Eq. [5] using the 
density of swimmers measured far away from the wall. 
In the limit of strong confinement (where R Ip), the 
bulk density tends to zero [s^, and thus the bulk swim 
pressure also vanishes. This results in the divergence of 
Eq. [¥]when Ip/R —>■ —oo. In this limit, nearly all particles 
are trapped at the wall, and the density distribution on 
the wall can be described in terms of the local curvature 



We now consider a sharp corner in a zig-zagging wall, 
with an opening angle 0 < 9 < tt, and examine the result¬ 
ing pressure near both the inner and outer corners. Near 
the inner corner, increased build-up of particles results in 
a locally higher pressure (as shown in Fig. [5^). In con¬ 
trast, close to the outer corner the local swim pressure 
is decreased. We note here that even for relatively long 
walls {L/lp ~ 20, there can still be a significant effect 
of the wall length on the measured force. We note that 
the total excess force on the wall (as compared to the 
force resulting from the global swim pressure Pq) is zero. 
This suggests that local irregularities in the geometry of 
a (globally flat) wall do not interfere with the observa¬ 
tion that the pressure of a gas of torque-free swimmers 
is well-defined . In we show for several angles 6 
the finite-size scaling for the inner corner for data ob¬ 
tained both from solving the Smoluchowski equation and 
from simulations. In Fig. [2j: we plot the total excess 
force Fex = Fact — LPq on the wall, extrapolated to the 
limit of small Ip, as a function of the opening angle. The 
behavior is well fitted by a function of the form 


that P{R) = Fo(l — Alp/R), with A ~ 0.836. Thus, 
the net pressure resulting from the forces on both sides 
of the object is given by Pnet{R) = 2AlpPo/R. If we 
assume (without loss of generality), that the curve ’w(s) is 
parametrized by length (i.e. |w'(s)| = I, then 1/R{s) = 
where cj){s) is the local angle of the curve with the 
x-axis, and Eq. [5] reduces to: 


Fact 


OAl P /^COS()-(si) -COS(/)(so)'\ 
P ° l^sin()!)(si) - sin())(so)y 


( 7 ) 


Thus, for curves with a smoothly changing, large radius 
of curvature, the total force on a curve is only dependent 
on the angles at the endpoints. For objects with sharp 
corners, the excess force from Fig. [5J: should be included 
in Fact as well. It is interesting to note that in both 
cases, the force on an object is independent of the size of 
the object. Assuming the persistence length of the active 
solvent is sufficiently small compared to the scale of the 
object, increasing the length of straight walls does not 
generate additional net forces. Scaling a curved object, 
on the other hand, changes the excess pressure on the 
wall due to the change in radius of curvature. However, 
in the limit where R ^ Ip this effect is cancelled out by 
the change in length of the curved region. 

In order to illustrate how these results can be used to 
predict the motion of obstacles in a bath of swimmers, we 
assume a simple setup where the object and swimmers 
are suspended in a three-dimensional fluid, but confined 
to two dimensions via e.g. optical forces. In this case, we 
can make use of slender body theory 38, ^ to determine 
the forces exerted on the object by the solvent. For exam¬ 
ple, the approximate drag force on a curved cylindrical 
object moving through a solvent moving with velocity v 
is given by: 


Fdrag — 




log(L/a) Jo 


I - iw'(s)w'(s) ) ds, (8) 


Fex(0) = 0.83Fo/pCot(6i/2). (5) 

The calculated forces can be used to predict the motion 
of large objects in contact with an active suspension. As 
an example, we consider a slender object with a contin¬ 
uously varying radius of curvature, for which the shape 
is parametrized by w(s). We assume that, compared to 
the persistence length Ip of the particles, the radius of 
curvature of the object R{s) is large everywhere, and its 
rate of change is small. The total force on this object is 
given by: 

F=f F(i?(s))n(s) |w'(s)| ds, (6) 

-'So 

where the prime indicates a derivative with respect to s, 
and n(s) is a vector that is locally normal to the surface 
of the object. In the limit of large radius of curvature 
R, the pressure on a curved wall is linear in 1/R, such 


where L is the total length of the curved object, a is the 
radius of the circular cross-section of the object, is the 
solvent viscosity, and I is a unit matrix. If the object is 
sufficiently large, we can assume that it will move slowly 
compared to the velocity of the active particles {v <^vo), 
and thus the forces exerted by the active particles will 
be independent of the object’s velocity. In this case, the 
velocity v can be calculated simply be setting Fdrag = 

Fact ■ 

In Fig. 131 we plot the predicted velocity for both 
a sharp V-shape and a continuously curved U-shape 
(parametrized by a cosine function) for a range of open¬ 
ing angles 9 and a fixed total length L. These shapes 
will be propelled forward, in the direction away from 
the “open” end of the shape. When the shape is ap¬ 
proximately a straight line {9 ~ tt), both shapes move 
approximately equally fast. However, there are differ¬ 
ences between the two velocities at smaller angles. This 
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FIG. 2. a) Local wall pressure near the inside and outside of a sharp corner of angle 6 = 7r/2, as a function of the distance 
X to the corner. The corners are part of a wall with a zig-zag geometry, as shown in the inset. Note that far away from the 
corner, we recover P — Pq. b) Net force near an inner corner of varying angle for different choices of the ratio of the wall 
length L and the persistence length Ip. The points are the results of Brownian dynamics simulations, the lines are obtained by 
numerically solving Eq. [3] c) Excess force on the wall next to a sharp corner of angle 0, in the limit where the lengths of the 
(straight) walls adjacent to the corner are much longer than the trajectory persistence length Ip. The blue and red points are 
obtained from the solution to Eq. [^and Brownian dynamics simulations respectively. The black line is the fit in Eq. [S] The 
force direction is always normal to the wall, and angles 9 > tv correspond to the outer corner. 



FIG. 3. Velocity u of a slender U-shape (red) or V-shape 
(blue) immersed in an active bath. The U-shapes are 
parametrized by a cosine function of varying amplitude (re¬ 
sulting in a varying angle 6 between the two ends), while the 
sharp shapes simply consist of two line segments at an angle 
9 (see pictures). The object length L is the same in all cases. 
The prefactor F = log(I//a)/47rpL, with fj. the viscosity of the 
solvent. The inset shows that the net force on each shape is 
approximately the same. 


framework can readily be extended to calculate the net 
velocity or angular velocity of an arbitrarily shaped slen¬ 
der object, assuming that all of its relevant length scales 
L ^ Ip. Interestingly, these differences purely result from 
differences in the friction the obstacle experiences from 
the solvent: as shown in the inset, the net force on each 
shape is the same for V- and U-shapes within our nu¬ 
merical accuracy. This can also be seen directly from the 
fact that the prefactor in Eq. [5] is approximately equal 
to A/2. 

Finally, we use the calculated forces on both curves 
and corners to predict “active depletion” forces: the in¬ 
teraction between two large disks immersed in an active 


solvent. In this case, we consider two disks of diameter 
CTc, immersed in a bath of particles with diameter cTa. The 
excluded volume for the active particles consists of two 
overlapping disks of diameter (ctc + cr(j)/2, with a center- 
to-center distance r. Where the two disks touch, the two 
sharp corners lead to a repulsive force between the parti¬ 
cles, which is counteracted by the active pressure on the 
rest of the dumbbell shape. The total depletion force is 
thus given by: 

Fdepi = 2E;x(6')sin(6l) - P{lp/ac)(7cSm{e/2), (9) 

where 9 = 2 arccos(r/crc) is the opening angle at the 
points where the two disks of excluded volume meet, and 
Pilp/cTc) and Fex{6) are given by Eqs. S] and El respec¬ 
tively. 

In Fig. Ell, we plot the active depletion force as a 
function of distance for several choices of the persistence 
length. In the limit of small Ip, the pressure on the 
curved walls dominates the forces on the corners, and we 
recover the two-dimensional Asakura-Oosawa depletion 
force, with an effective temperature ksTes = C,v^l2Dp 
(such that Pq is equivalent to an ideal gas pressure). We 
obtain good agreement between simulations and the pre¬ 
diction of Eq. Elfor small values of Ipjac- We can obtain 
an effective interaction potential, Uef{{r) by integrating 
the obtained forces: 

Ueff(r) = - f dr'Fact(?'')- (10) 

J (Jc+iya 

where we have chosen the interaction potential to be zero 
at r = (Tc + CTa. We plot the effective potentials in Fig. 01. 
Note that despite the divergence in the force (Fig. Eli), 
the interaction energy remains finite, indicating that the 
colloids can in principle approach each other. 

Note that for distances larger than r = cFc + era, the 
active particles can suddenly pass between the two disks. 
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FIG. 4. a) Active depletion force between two disks of diameter immersed in an active bath (particle diameter aa = O.Ictc) 
for various values of the persistence length Ip. The lines result from numerical calculations, the points from simulations. The 
black dashed line is the depletion force in the two-dimensional Asakura-Oosawa model, with the temperature set to the effective 
temperature fcsTeff = (^Vol2Dr. Note that a singularity occurs when r = ac + (Ja. b) Schematic representation of the disks 
(red), depletion zones (white), and the relevant length scales, c) Effective interaction potential C/eff corresponding to the forces 
plotted in panel (a). 


leading to a discontinuous change in the force. Although 
we expect the depletion force to decay rapidly to zero 
beyond this distance, we have not closely investigated 
this regime, as it inherently involves geometrical features 
with a length scale smaller than Ip. However, simulations 
indeed show a steep drop to zero in the force at slightly 
larger distances. Moreover, in a recent study on active 
depletion forces [2l|, simulations for relatively weak ac¬ 
tive propulsion (i.e. small Ip) indeed revealed a peak in 
the depletion force near r = CTc + <7o, and a decay to zero 
at larger distances. Note that more complex effects, such 
as geometric shielding, can occur when the persistence 
length of the particles is large compared to the separa¬ 
tion of the two obstacles. In particular, simulations of 
swimmers confined between two walls [l9| showed oscil¬ 
latory behavior of the interaction force in this regime. 


In conclusion, we have generalized the concept of swim 
pressure to non-flat walls, and calculated the curvature- 
and corner-induced corrections to the bulk swim pressure 
both numerically and via particle-resolved simulations. 
This approach yields a coherent framework to obtain the 
forces and torques on passive particle of arbitrary shape 
and to derive a simple approximation for the activity- 
induced depletion interaction. 


Our approach can in principle be generalized to three 
spatial dimensions as well as to incorporate more de¬ 
tails of the hydrodynamic swimmer-wall interactions. It 
would further be interesting to study the motion of pas¬ 
sive carriers in a density gradient of active particles. 
Moreover, an important extension of this work would 
be to consider an active background outside the dilute 
regime, where the interactions between the active parti¬ 
cles can lead to clustering, ordering, and phase separation 

i3-Ii5|. 


APPENDIX: COMPUTATIONAL DETAILS 


In order to numerically obtain the density distribution 
of active particles near an obstacle, we find the steady- 
state solution of the the Smoluchowski equation for the 
orientation-dependent density distribution p(x, 


dp{^, <(>, t) 
dt 


= ■ Vp(x, (l),t) +D, 


9(^2 


= 0 . 

( 11 ) 

Here, x indicates the particle position, (j) denotes the 
particle orientation, and v((^) = rip (cos (^, sin <()) is the 
self-propulsion velocity of a particle with orientation </>, 
where vq is the swim speed. Dp is the orientational diffu¬ 
sion coefficient. To find the steady-state solution to this 
equation, we divide the system area up into Uxy cells, and 
divide the set of possible orientations into equal inter¬ 
vals. We assume that the cells are sufficiently small, such 
that in each cell p(x, (p) is approximately constant over 
the area and orientation range of the cell, and calculate 
the total flow of particles into and out of each cell. For 
each cell, this particle current can be written as a linear 
combination of the densities of the cell and its neighbors. 
Setting the rate of density change for each cell to zero 
yields a complete set of linear equations for the density 
Pi^k for each cell i and orientation k, of the form: 


= vo^Tij{k)pj{(j)k) + Dr^Rklii)p^i^l) = 0 , 

j=i 1=1 

( 12 ) 

where the translation matrix Tij{k) is the Uxy x Uxy ma¬ 
trix representing the flow of particles with orientation k 
between cells due to their active motion. Similarly, the 
rotation matrix Rki{i) represents the rotational diffusion 
in each cell i. In order to obtain the elements of the trans¬ 
lation matrix, we determine the particle current Jc,e out 
of cell c through edge e with edge length and norm ne 
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using an upwind scheme, yielding: 


v{cj}k) ■ nelePc,k a v{cj}k) ■ Ue > 0 

v{(j)k) ■ ne^ePneighbor.fc Otherwise 


(13) 


where 4>k = 2kTr/n^, and Pneighbor.fc indicates the density 
of the cell neighboring cell c through edge e. Combining 
all particle currents yields the translation matrix: 


/ ^Pi,k \ 

V dt 


E 

e€edges(z) 



^0 E 

j^cells 


(14) 

(15) 


where Ai is the surface area of cell i. The rotation ma¬ 
trix Rki{i) is obtained by numerically approximating the 
second-order derivative with respect to (f) in Eq. fTTl 


dPi,i 

dt 


_ jj [Pi,k+1 Pi,k] [Pi,k Pi,k—l] 


{27Tfn^)^ 

— Rr ^ ^ Rkl i'^)Pi,l 


(17) 


Using Eq. the translation and rotation matrices 

define a full set of linear equations for the density distri¬ 
bution in the system, which can be solved numerically. 
The hard walls in the system are represented by zero-flux 
boundary conditions on parts of the mesh edges. For the 
other mesh boundaries, we make use of either reflecting or 
shifted boundary conditions, depending on the symmetry 
of the geometry. For the case of curved walls, we make 
use of the rotational symmetry that the density profile 
exhibits both on the inside and the outside of a circular 
wall. In this case, only one row of cells is needed along 
the radial direction, and translation along the azimuthal 
direction can be considered as a change in orientation in¬ 
stead of a change in position. In the case of walls with 
sharp corners, we implement a wall with a zigzagging ge¬ 
ometry, where we make use of mirror symmetry in order 
to reduce the size of the mesh. Additionally, to reduce 
numerical error, we make use of the knowledge that the 
total wall pressure should equal the bulk swim pressure 
(Eq. El). 
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